Magnetic interference pattern in planar SNS Josephson junctions 
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We study the Josephson current through a balhstic normal metal layer of thickness D on which two supercon- 
ducting electrodes are deposited within a distance L of each other. In the presence of an (in-layer) magnetic field 
we find that the oscillations of the critical current Ic{^) with the magnetic flux $ are significantly different from 
an ordinary magnetic interference pattern. Depending on the ratio L/D and temperature, /c(<l?)-oscillations can 
have a period smaller than flux quantum $o, nonzero minima and damping rate much smaller than 1/$. Similar 
anomalous magnetic interference pattern was recently observed experimentally. 

PACS numbers: 74.78.FK, 74.50.+r 



I. INTRODUCTION 

Existence of a supercurrent in a Josephson junction is the 
manifestation of the interference between the macroscopic 
wave functions (superconducting order parameters) of the two 
contacted superconductors. The quantum interference can be 
modulated by an external magnetic field applied to the junc- 
tion. As the result the critical(maximum) supercurrent Ic{^) 
shows a well-known Fraunhofer diffraction pattern-like de- 
pendence on the magnetic flux $ penetrating the junction area. 
In a superconductor-insulator-superconductor (SIS) junction 
the critical current. 



sin(7r$/$o) 
7r<i>/$o ' 
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oscillates with the period of flux quantum $o = Trfic/e and 
an amplitude decreasing as !/<!> The main features of the 
effect, i.e. damped oscillations of Ic with the magnetic flux, 
take place in other types of Josephson weak links, however 
the detailed behavior, including the period of the oscillations 
and the rate of damping, depends on the geometry as well as 
the nature of the weak link. 

In a wide SNS (N being a normal metal layer) junc- 
tion Ic has a similar magnetic interference pattern as SIS 
systems^. On the other hand, Heida et at.,— investigating S- 
two-dimensional-electron-gas-S (S2DEGS) junctions of com- 
parable width W and length L, have measured a 2$o period- 
icity of the critical current instead of the standard <I>o period- 
icity. The first explanation of this finding was due to Barzykin 
and Zagoskin^, who considered a S2DEGS junction with per- 
fect Andreev reflections at NS interfaces and both absorbing 
and reflecting lateral boundaries, obtained a 2<I>o periodicity 
in the limit L/W oo (i.e. in the limit of the point-contact 
geometry). Later, Ledermann et al.,^ considered more real- 
istic reflecting boundaries at the edges and found that in the 
limit of strip geometry (L/W ~ 1) the periodicity of the crit- 
ical current changes from <f>o to 2<I>o as the flux through the 
junction increases. In general, increase of the periodicity was 
attributed to the nonlocality of the supercurrent density in hy- 
brid NS structures. 

In this paper we report on a new type of magnetic interfer- 
ence pattern in a planar SNS junction. The SNS junction stud- 
ied below consists of a thin normal metal layer of thickness D 



on which two superconducting electrodes are deposited in a 
distance L of each other (Fig.[T]|. An external magnetic field 
is applied in plane of N-Iayer perpendicular to the direction 
of the Josephson current flow. Such a Josephson setup was 
recently used in the experiment by Keizer et ali^ to inves- 
tigate the Josephson supercurrent through a S-half-metallic- 
ferromagnet-S (SHMFS) junction. They used NbTiN super- 
conducting electrods on top of a thin layer of Cr02, which is 
a fully polarized (half-metallic) ferromagnet, i.e. it supports 
only one spin direction of electrons. Surprisingly, Josephson 
supercurrent was detected for the junction length L ^ 300nm 
in spite of strong pair breaking in Cr02 which is expected to 
suppress all singlet superconducting correlations. Ref. 7 also 
reports measurement of the magnetic interference pattern with 
magnetic field applied in plane of HMF-layer. It is found that 
Ic{^) oscillates with the in-plane flux <i> with a period of or- 
der $0. In contrast to a standard magnetic interference pattern 
dl.lt , Ic has nonzero values at the minima and the amplitude 
of the oscillations decreases rather slowly compared to l/(i>. 

The long range superconducting proximity in HMF can be 
explained in terms of the triplet superconducting correlations 
generated at spin-active HMFS-interfaces as the result of in- 
terplay between the singlet superconducting correlations and a 
noncollinear magnetization inhomogeneity*^-^. While the sin- 
glet superconducting correlations are destroyed over a short 
distance of the Fermi wave-length, the triplet components 
can survive over distance of the order of the normal coher- 
ent length — wf/ZttT, with v-p being the Fermi velocity. 
In Ref. 7 by investigating the critical supercurrent for differ- 
ent distances between the electrodes it was concluded that the 
length dependence is the same as in nonmagnetic SNS junc- 
tions. This strongly suggests indeed that triplet correlations 
are responsible for the observed Josepshson current. 

The aim of the present paper is to investigate /c($) in such 
a planar SNS junction. Whereas this problem is interesting 
by itself, we also believe that it is relevant for the understand- 
ing of the experimental results of Ref. 7. Indeed, penetration 
of triplet correlations in SHMFS junction is similar to an or- 
dinary (singlet) superconducting proximity in SNS systems, 
i.e. both decay exponentially within the length scale ^at. By 
noting this fact we will use the quasiclassical Green's func- 
tion formalism to investigate the magnetic flux dependence 
of the supercurrent in the corresponding planar SNS junction 
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FIG. 1: (Color online). Schematic of the SNS junction. The current 
flow between two superconductors (81,2) through a Normal layer 
(N). The NS interfaces are perfectly transparent. 

(see Fig. [TJ. We find that the magnetic interference pattern 
is significantly different from that of the standard one. The 
period of the oscillations can be smaller than $0 depending 
on the length-to-thickness ratio L/D. The period tends to $0 
at higher magnetic fluxes and also for very large L /D. We 
also obtain the two anomalous features observed in the ex- 
periment: the amplitude of the oscillations has a rather slow 
decrease with $ compared to the standard SIS case (II. lb . and 
the critical current as the function of the flux, /c($), at low 
temperatures can have finite minima when the total flux $ is 
integer or non-integer multiplies of the $0 depending on the 
period of the oscillations. 

In Section we introduce our model of a ballistic planar 
SNS contact and present solutions of the Eilenberger equation 
for the quasiclassical Green's functions of a given electronic 
trajectory. Introducing the effect of the in-plane magnetic field 
through a gauge invariant phase we obtain the expression of 
the critical supercurrent as a function of the magnetic flux. 
Section Ulllis devoted to the analysis of the Ic{^) in terms of 
L/D for different temperatures. In Section|IV]we present the 
conclusion. 



II. JOSEPHSON CURRENT IN SNS JUNCTION WITH AN 
IN PLANE MAGNETIC FIELD 

In this section, we calculate the Josephson cuiTent for a 
clean SNS junction in the presence of an external magnetic 
field H. The setup is schematically shown in Fig. [T] It con- 
sists of a normal metallic (N) layer of the thickness D and 
width W on which two superconducting electrodes are de- 
posited in a distance of L. This planar Josephson structure 
was studied experimentally in Ref. 7 with a half metallic 
N layer. A phase difference between order parameters of 
the superconductors drives a Josephson supercurrent through 
the parts of N layer underneath the superconductors and the 
junction N part {—L/2 < y < L/2). The magnetic field 
H = —Hx is applied in the plane of N layer perpendicu- 
lar to the direction of supercurrent flow. We consider a clean 
structure with all dimensions L,W,D being smaller than the 
electronic impurity mean free path ^imp and ideally transpar- 



ent NS-interfaces. At the same time, the Fermi wavelength 
is small compared to L,W,D, and the superconducting co- 
herence length ^0 = i'F/2Ao(r = 0) (we use the system of 
units with h= k-^ ~ 1). Under these conditions the electronic 
properties of the system can be derived from the Eilenberger 
equations for the semiclassical matrix Green's function g, 

-i;Fn.Vt = c^„[f3,^] + [A(r),^], (2.1) 
The matrix Green's function 

-(/!:-!)■ 

where the normal g and anomalous / Green's functions de- 
pend on the Matsubara frequency cj„ = 7rT(2n + 1), on the 
coordinate r, and the direction of motion n; ti (i=l,2,3) de- 
notes the Pauli matrices in the Nambu space, and the matrix 

A(r) - ( « ^ 

represents the superconducting order parameter A(r). The 
matrix Green's function g^^ satisfies the normalization con- 
dition, 

t' = l, ^L+/--./l=l' (2-2) 

where /J^^(r, n) ~ /* (r, — n). In components, the Eilen- 
berger equations have the form 

-v^n.^g^^^ - A(r)/t^-A*(r)^„, (2.3) 
-v^n.yU^ = 2c^„/^„ -2A(r)g^^. (2.4) 

We solve these equations along an electronic quasiclassical 
trajectory (shown in Fig. [TJ which is parameterized by — cx) < 
T < oc"^'". Assuming a weak external magnetic field we 
neglect its effect on the orbital motion of the quasiparticles. 
The magnetic field then will have only phase effect which we 
will include by introducing the gauge invariant phase as 

27r /"^ 

= V'o - -;- / A.ndr, (2.5) 
'^0 Jo 

where ipo is the phase difference between two superconduc- 
tors in the absence of the magnetic field, the second part is the 
phase accumulated by the quasiparticle on the trajectory due 
to the magnetic field H = V x A, with £ being the length of 
the trajectory inside N-layer The vector potential is taken to 
be A = —Hzy. 

A typical trajectory consists of three parts: the part ex- 
tended from bulk of Si (r = —00) to a point at NSi -interface 
T — 0, the part inside N layer (0 < t < i) which extends from 
a point at NSi-interface to a point at NS2-interface and the 
last part which extends from the point r = ^ to the bulk of S2 
(r = 00). Our equations are supplemented by the boundary 
conditions which determine the values of Green's functions in 
the bulk Si and S2, 

f Ao(T) exp(iiy.i^2) 

Juj„[T = Too) = , (2.6) 

g^Jr^T^) = ^, (2.7) 
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where r2„ = -^/Aq + ujf-^, and Ao(T) is the temperature- 
dependent superconducting gap. We neglect the variation of 
the order parameter close to NS-interfaces inside the super- 
conductors and approximate the order parameter by the step- 
function, A(r) = Aoe"^ie'(r) + Aoe*'^^6'(T - L). We can 
obtain the Green's function g^^ on a trajectory inside N that 
is constant. It depends only on the length of that trajectory £ 
and the phase difference ipo, which is given by 

g^^ = ta.nh[uJn£ / vp + i(po/2 + arcsinh (aj„/Ao)] . (2.8) 

Note that in the presence of the magnetic field, the phase dif- 
ference (fio is replaced by if (see Eq. ( 12.5b ). 

The supercurrent density can then be obtained by averag- 
ing ( I2.8I 1 over all different possible classical trajectories. This 
corresponds to an averaging over Fermi velocity directions. In 
the presence of the planar magnetic field, we find 

j(r) = eA/'(0)TwF^ J nlm g^Ji, vp) sin 9 dedc/) . (2.9) 



For the first term of the right side, we can write 



A.ndr — —H / zdy . 
Jo 



Since the equation for this segment of the trajectory is z 
D{1 — y{D ta.n 9 sin (f))^^ , we get for the integral 



A.ndr 



HDyo 



,yo 



2D tan 9 sin < 



(2.13) 



where yo is shown in Fig. [T] Similarly, we find identical re- 
sults for the integrals of the vector potential over the other 
segments. Therefore for a trajectory with the length £, the 
phase induced by the planar magnetic field proportional to 



Here A/^(0) is the density of states at the Fermi surface and 
Im g^^ {£, vp) denotes the imaginary part of the normal com- 
ponent of the matrix Green's function, given by 



Al{T) sin LP 



(fll+Lol) cosh X + 20„ w„ sinh x + Aq cos (yS ' 



X 



(2.10) 
(2.11) 



A.nrfr = 



-HDNyo 



(2.14) 



where N is the number of the triangles (each triangle con- 
sists of two segments) for the trajectory of the length £ passing 
through the point z, N ^ [L/{4Dta.n9sin<j)) - z/{2D)] + 
[L/(4£)tan6'sin0) + z/{2D)] + 2, with the square brackets 
denoting the integer part. 

Thus, for the phase difference, we obtain 



To calculate the integral of the vector potential along the 
quasiclassical trajectories, we split it into the segments as 
shown in Fig. [1] 



^0 



$ 2iVi:)tan6'sin( 
L 



(2.15) 



A.ndr = / A.ndr + / A.ndr + ... 

"'0 Jti 

r„-i pi 

A.ndr + / A.ndr . (2.12) 

I 



Here $ = HDL is the total flux through the junction. Sub- 
stituting this into Eq. ( 12.91 ) and taking the y-component of the 
cmTent, we obtain the final expression 



T 



E 



D /.TT 



A2(T)sin +7r^^(l -a;2)i/2sin</.) {I ~ x^^/^dx sin (l)d(l)dz 
Jo "/-i (f]2 _|_^2)coshx + 2f]„tj„sinhx + A2cos(^(po + 7r^f(l -a::2)i/2sin(/)y 
I{^o) 



0<lpq<2ti 



lo 



(2.16) 



with the notations cos 6* = x and Iq — 2evpJ\f{Q)TcW^o, 
and W being the width of the normal layer in the a;-direction, 
which for a wide junction is taken to be much larger than D 
and L. The length £ of the quasi-classical trajectory equals £ — 
2ND/x. For low temperatures, T ^ vp / L, the summation 
over Matsubara frequencies in Eq. (12.16b can be replaced by 



the integration. 



2ttT 



]{...) ^ j{...)duj-^ y"(...)Ao(r) cosh^d/i 



where = Aq sinh ji. For a long SNS junction (L ^ ^0) and 
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FIG. 2: The critical current dependence on the external magnetic 
field, applied in the plane of the normal metallic layer for different 
temperatures, L/D ^ 100 and D/^o = 1/30 . 



low temperature, Eq. ( |2.16l l is simplified as 

l-D f-TT fl 



lo 



^0 



2D f^o 



^0 



£2 



J sin(m(p) 



, (2.17) 



rn—l 



where E{(p^ i) is the Fourier Series, given by (see Ref. 

.(2.18) 



27r 2 



III. DISCUSSION AND RESULTS 



Equation ( 12.161 ) expresses the magnetic interference pattern 
/c(<i>) of a ballistic SNS junction in the presence of an in-plane 
magnetic flux $. In this section we analyze Id^) in terms of 
the length-to-thickness ratio L /D and the temperature T for 
D/^o = 1/30. 

Let us start with analyzing the case of very large L/D. Fig- 
ure |2]a,b shows oscillations of /c(<&) for L/D = 100 and at 
low (T = O.lTc) and high (T = O.QSTJ temperatures, re- 
spectively. At low temperatures the the critical current goes 
through nonzero minima at finite fluxes. The amplitude of 
supercurrent minima decrease with $ and drops to zero at 
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FIG. 3: The same as Fig.[2]but for L/D = 10. 



$ ^ $0- Compared to an ordinary magnetic interference pat- 
tern, the oscillations are weakly damped since their amplitude 
decreases with $ much slower than as 1/$. With increasing 
temperature, the amplitude of the oscillations decreases. Also, 
the minimal values of the supercurrent decreases and vanishes 
as T Tc, where Ao(T) < T. Note that at both low 
and high temperatures the period of oscillations varies from 
0.92$o (first minima) at low magnetic fluxes to $o at high 
fluxes. The result that the period of oscillations is tempera- 
ture independent comes from the fact that the gauge invariant 
phase in the argument of the sine and cosine functions in Eq. 
(I2.I6I 1 does not contain any temperature-dependent factors. 

Figure |3] presents the magnetic interference pattern for a 
lower L/D —10 and at the same temperatures as Figure |2] 
From these plots we see that lowering L /D has two main ef- 
fects. First, the rate at which the amplitude of Ic oscillation 
decreases with $ increases. Second, the period of oscillations 
at both low and high temperatures becomes smaller ~ 0.72$o 
at small fluxes. Increasing $ the period increases up to 2± $o- 
Again as in Figs. |2]the value of supercurrent at the minima 
vanishes as the temperature approaches Tc- 

Still lower period of oscillations at small fluxes can be 
reached at low values of L/D. This is illustrated in Fig. 
|4] where the magnetic interference pattern is presented for 
L/D = 5. Clearly the decay is close to the ordinary pat- 
tern ( II. Il l, i.e. 1 /$, and the period can be as small as half the 
flux quantum. 

The existence of the nonzero minima in the oscillations of 
/c($) is related to the non-sinusoidal phase dependence of 
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FIG. 4: The same as Fig.|2]but for L/D = 5. 



the Josephson current ( 12.161 ) which is more pronounced at low 
temperatures T <^ Tc- We have found that Ic{^) undergoes 
a change of sign at a nonzero minimum. At such point the 
amplitude of the first harmonic (oc sin ip^) of the Josephson 
current vanishes and Ic{^) is determined by the amplitude 
of the higher (mainly second) harmonics which change sign 
upon crossing the minimum. Similar effect was found before 
in ferromagnetic Josephson junctions (see Ref. ITslflSl and ref- 
erences therein). As T approaches Tc the ratio Aq (T)/T goes 
to zero and the current-phase relation (Eq. ( |2.16t ) becomes si- 
nusoidal and consequently the nonzero minima disappear. 

The dependence of the period of the oscillations on the 
magnetic flux and the geometry can be understood in terms 
of the difference between the magnetic flux $£ (Eq. (|2.14t ) 
enclosed by a trajectory of length t, and half of the flux $ 
penetrating through the area DL. The difference comes from 
the fact that a trajectory which does not pass through the edges 
of S contacts has extra parts in N-layer which lie outside the 
area (See Fig. [D. Writing = $/2+(5<l5^, the difference 
5^1 vanishes only for the trajectories which pass through the 
edges of S contacts. A finite averaged ((5$^ ) over different tra- 



jectories means that the period of /c(<I>)-oscillations, obtained 
from Eq. ( |2.16t . differs from $o- We note that in the limit 
of thin N-layer L ^ D or high magnetic fluxes $ $o the 
contribution of the trajectories which are not passing through 
the edges is negligibly small in the interference structure and 
the period of the oscillation approaches $0 (See Figs. |2]|4|. 

In contrast to L 3> -D case, for smaller L/D the contribu- 
tion of the trajectories not passing through the edges is impor- 
tant. Because of having larger length, a trajectory not passing 
through the edges has bigger contribution in gaining the effect 
of the magnetic flux as compared to the corresponding trajec- 
tory (having the same orientation 6 and (p) passing through the 
edges. Therefore we expect that the effect of magnetic flux is 
more pronounced for thicker N layers compare to the thinner 
ones, which explains why the decrease of the amplitude of 
/c(<i>) -oscillations with <1> is faster for smaller L/D. 



IV. CONCLUSIONS 

In conclusion, we have studied Josephson effect in a SNS 
structure made of a thin ballistic N-layer of thickness D on 
which two superconducting electrodes are deposited at the 
distance L between each other. A magnetic field is applied in 
plane of the N-layer which modulates the superconducting in- 
terference and leads to a decaying oscillatory variation of the 
critical supercurrent /c($) with the magnetic flux $. Using 
the quasiclassical Green's functions approach, we have shown 
that such a magnetic interference pattern has three main dif- 
ferences with that of an ordinary pattern. First, at low temper- 
atures the oscillations of the critical current /c(<&) go through 
the minima at which the supercurrent has nonzero values. Sec- 
ond, for a large L/D the amplitude of the quasi-periodic os- 
cillations of Ic{^) decays at a rate which is much slower than 
1/$. Third, at low magnetic fluxes the oscillations can have a 
period smaller than the magnetic quantum flux $0 depending 
on L/D. These features have been experimentally observed 
recently^. 
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